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Abstract. We report the results of a search for a point X-ray source (stellar remnant) in the southwest protrusion 
of the supernova remnant G 315.4—2.30 (MSH 14— 65, RCW86) using the archival data of the Chandra X-ray 
Observatory. The search was motivated by a hypothesis that G 315.4—2.30 is the result of an off-centered cavity 
supernova explosion of a moving massive star, which ended its evolution just near the edge of the main-sequence 
wind-driven bubble. This hypothesis implies that the southwest protrusion in G 315.4—2.30 is the remainder 
of a pre-existing bow shock-like structure created by the interaction of the supernova progenitor's wind with 
the interstellar medium and that the actual location of the supernova blast center is near the center of this 
hemispherical structure. We have discovered two point X-ray sources in the "proper" place. One of the sources 
has an optical counterpart with the photographic magnitude 13.38 ± 0.40, while the spectrum of the source can 
be fitted with an optically thin plasma model. We interpret this source as a foreground active star of late spectral 
type. The second source has no optical counterpart to a limiting magnitude ~ 21. The spectrum of this source 
can be fitted almost equally well with several simple models (power law: photon index = 1.87; two-temperature 
blackbody: kTi =0.11 keV, Ri = 2.34 km and kT2 = 0.71 keV, R2 = 0.06 km; blackbody plus power law: 
kT = 0.07 keV, photon index = 2.3). We interpret this source as a candidate stellar remnant (neutron star), while 
the photon index and non-thermal luminosity of the source (almost the same as those of the Vela pulsar and 
the recently discovered pulsar PSR J 0205-1-6449 in the supernova remnant 3C 58) suggest that it can be a young 
"ordinary" pulsar. 

Key words. Stars: neutron - ISM: bubbles - ISM: individual objects: G 315.4-2.30 (MSH 14-65, ROW 86) - ISM: 
supernova remnants - X-ray: stars 



1. Introduction 

G 315.4-2.30 (MSH 14-65, RCW86) is a bright (in both 
radio and X-ray) shell-like supernova remnant (SNR) with 
a pec uliar protrus ion to the southwest (e.g. Dickel et al. 
2001 , Vink et al. |l997 ). This protrusion encompa sses a 
bright hemi spherical optical nebula (Rodger et al. 1960 , 
Smith 1997 see also Fig. ^). The characteristic angular 
size of the SNR is about 40 , that at a distance to the 
remnant of 2.8 kpc (Rosado et al. 1996|) corresponds to 
~ 32 pc. The radius of the optical nebula is ~ 2 (or ~ 1.6 
pc). A collection of observational data points to the young 
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age ( few thousand years) of the SNR (see e.g. Dickel et al. 
^OOll) . 



Vink et al. ( |1997| ) put forward the idea that the SNR 
G 315.4—2.30 is the result of a supernova (SN) explosion 
inside a pre-existing wind-driven cavity (see also Dickel 
et al. |2001| , Vink et al. |2002D and noted that the elon- 
gated shape of the SNR resembles that of a wind-driven 
cavity created by a moving massive star (see Weaver et 
al. |1977| and Brighenti & D'Ercole |1994| for details). On 



the other hand, it is believed that the origin of the south- 
west protrusion is due to the interaction of the SN blast 
wave with a density enhancement in the ambient interstel- 
lar medium. However, as was correctly noted by Dickel et 
al. (2001), the density enhancement (e.g. a high-density 
cloud) should result in a concave dent in the shell, not in 
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a protrusion. Dickel et al. ( 2001 ) also suggested that this 
protrusion "is perhaps the key to what is going on" . 

We agree with the idea that G 315.4—2.30 is a diffuse 
remnant of an off-centered cavity SN explosion and supple- 
ment this idea by the suggestion that the massive SN pro- 
genitor star exploded near the edge of the main-sequence 
bubble. This suggestion implies that the southwest pro- 
trusion is the remains of a bow shock-like structure cre- 
ated in the interstellar medium by the post- main-sequence 
winds (see Sect. 4 and Gvaramadze 2002; cf. Wang et al. 
1993| ) and that the SN exploded near the center of this 



hemispherical structure. Given the youth of the SNR and 
assuming a reasonable kick velocity for the stellar rem- 
nant, one can expect that the stellar remnant should still 
be within the protrusion. Motivated by these arguments, 
Gvaramadze ( ^002 ) searched for a possible compact X- 
ray source to the southwest of G 315.4—2.30 using the 
i? 05^ T" archival data, but the moderate spatial resolution 
of the ROSAT PSPC precluded detection of point sources 
against the bright background emission of the SNR's shell. 

In this paper we report the discovery of two 
point X-ray sources near the center of the hemispher- 
ical optical nebula using the archival Advanced CCD 
Imaging Spectrometer (ACIS) data of the Chandra X-Ray 
Observatory. We interpret one of the sources as a fore- 
ground active star of late spectral type, and the second 
one as a candidate stellar remnant (neutron star). 

Note that Vink et al. ( ^000 ) also reported the dis- 
covery of a point X-ray source in the southwest half of 
G 315.4—2.30, at about 7 from the geometrical center of 
the SNR. We recall that Vink et al. ( |1997D suggested that 
the SN blast wave in this SNR takes on the shape of the 
pre-existing elongated cavity (created by the stellar wind 
of the moving SN progenitor star). However, the initially 
spherical shape of the wind-driven cavity could be signif- 
icantly affected by the stellar motion only if the massive 
star reaches the edge of the cavity and the stellar wind 
starts to interact directly with the ambient interstellar 



medium (e.g. Weaver et al. 1977); this implies that the 
SN explodes near the edge of the future (young) SNR (see 
Sect. 4 ). However, the source discovered by Vink et al. 
(|2000D is located too far from the edge of G 315.4-2.30. 
Moreover, the spectral characteristics of the source and 
the presence of a possible optical counterpart suggest that 
it is an active star rather than the stellar remnant associ- 
ated with the SNR (Vink et al. |2000[ see also Sect. 3.1). 



2. Point X-ray sources 

2. 1 . X-ray images 

Fig. |l| shows the Chandra (0.7—2 keV) image of the highly 
structured southwest corner of the SNR G 315.4—2.30. 
Two point X-ray sources are clearly visible at the north- 
west of the image at a2ooo — 14'^40™31!33, (^2000 ~ 
-62°38'22'.'8 and agnnn = 14i^40'"31!05, (5, 



Fig. 1. Chandra image of the southwest region of the SNR 
G 315.4—2.30. Both X-ray sources are surrounded by a 
circle. North is up and east is to the left. 

Fig. 2. The enlarged view of the region around the point 
X-ray sources, labelled S (southern) and N (northern). 

Fig. 3. The DSS-2 image of the optical nebula southwest 
of the SNR G 315.4-2.30. The position of both point X- 
ray sources is surrounded by a circle. 

Fig. 4. The enlarged DSS-2 image of the region around 
the point X-ray sources. The circles are centered at the 
position of X-ray sources. 



function at the off- axis angle 2'.5. Fig. ^ shows a close- 
up of the region around the sources, labelled S (southern 
source) and N (northern source). 

2.2. Optical images 

Fig. H shows the image of the hemispherical optical nebula 
to the southwest of the SNR G 315.4-2.30 from the Digital 
Sky Survey (DSS-2, red plates). Fig. ^ shows an enlarged 
view of the region around the point X-ray sources. The 
source S is in good positional agreement with a point- 
like object immersed in a diffuse optical filament. This 
object is indicated in the HST Guide Star Catalog as non- 
stellar, with a photographic magnitude of 13.38 ± 0.40. 
We believe, however, that it was misclassified due to the 
effect of the background diffuse emission and in fact it is a 
star. The source N has no optical counterpart to a limiting 
magnitude ~ 21. 

2.3. Spectral analysis 

The spectra of the point sources were extracted from cir- 
cular regions with radii of 1.5 arcsec. The background 
spectrum was taken from a circle with a radius of 4.5 
arcsec at 6.5 arcsec northeast of the point sources. The 
PSF model shows that the fraction of the source flux in 
the background region is negligible. The spectral modeling 
was performed in the 0.5—10 keV energy range using the 
XSPEC spectra fitting package. The estimates of the in- 
terstellar absorption, iVn, towards G 315.4—2.30 are quite 



uncertain (cf. e.g. Vink et al. 1997 with Vink et al. 2002), 
so the spectra were fitted with A^h as a free parameter. 
In one case, however, we used the fixed value of (see 
below) . 

The spectrum of the source S is shown in Fig. H; the 
solid line represents the best fit optically thin plasma 
model, with a temperature = 0.69 ± 0.05 keV, abundance 
-62°38'l6'.'9. For both sources we measured a FWHM of = 0.07 ± 0.02, A^h = (1-2 ± 0.3) x lO^^ cm~^ luminosity 
I'.'l, which is consistent with the Chandra point spread (in the 0.5-10 keV energy range) ~ 2.8 x lO'^^ ergs~^, and 
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1 2 5 

channel energy (keV) 

Fig. 5. The background-subtracted Chandra ACIS spectrum from the source S. The soHd line represents the best fit 
opticaUy thin plasma model. 



X^/d.o.f. =21.6/25. The blackbody or power law models 
give unacceptable fits. 

The spectrum of the source N (shown in Fig. ^ can 
be fitted almost equally well with a power law (PL), 
two-temperature blackbody (BB-f-BB), or blackbody plus 
power law (BB-I-PL) models. A simple BB model gives 
unacceptable fits. The best fit spectrum predicted by the 
BB-I-PL model is shown in Fig. | by the solid line. 

We note that the best fit PL model requires A^h = 
(we consider this fact as an indication of the galactic 
origin of the source N), although models with Ah up to 
1.5 X 10^^ cm"'^ (presented in Table 1) also give acceptable 
fits. 

The results of a spectral analysis for the source N are 
summarized in Table ^ 

3. Discussion 

3.1. Source S 

The soft X-ray spectrum of the source S and the existence 
of the candidate optical counterpart strongly suggest that 
this source is an active star. Assuming that S is indeed a 
star and given the observed X-ray to optical fiux ratio, one 
can exclude the possibility that this star is of the OB-type. 



Therefore it cannot be a member of a group of OB stars 
located in the direction of G 315.4— 2.30 n early at the same 
distance as the SNR (Westerlund [19661) . Most likely the 
source S is a foreground active star of late spectral type. 
A follow-up study of this object could clarify its nature. 

Note that the point X-ray source discovered by Vink 



et al. (2000) has much in common with the source S. The 
spectra of both sources are soft and could be fitted with 
the optically thin plasma model with depleted abundances 
(typical of late-type stars; e.g. Giampapa et al. 1996). 
Both sources have candidate optical counterparts with al- 
most the same photographic magnitudes. In both cases 
the X-ray to optical fiux ratios suggest that the optical 
objects are late-type stars. Evidence for long-term vari- 



ability of the former source (Vink et al. 2000 ) also suggest 
that it is an active star. 



3.2. Source N 

The spectral characteristics of the source N coupled with 
the absence of an optical counterpart allow us to consider 
it as a candidate stellar remnant. We now discuss this 
possibility. 
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channel energy (keV) 

Fig. 6. The background-subtracted Chandra ACIS spectrum from the source N. The sohd hue corresponds to the best 
fit blackbody plus power law model. 

Table 1. Spectral fits to the point X-ray source N. 



Model 


Photon index 


Temperature, 
keV 


10^^ cm~^ 


Luminosity, 
10^^ ergs"^ 


/d.o.f.- 


PL 


1.87 ±0.09 




1.5 (fixed) 


0.44 (0.5 - lOkeV) 


51.5/27 


B+B 




0.11 ±0.03 
0.71 ±0.05 


4.3 ± 1.8 


1.32 (bolometric) 
0.35 (bolometric) 


47.3/25 


BB+PL 


2.31 ±0.30 


0.070 ± 0.003 


9.2 ±3.4 


112 (bolometric) 
0.67 (0.5 - lOkeV) 


38.9/25 



Note that a single point near the Chandra mirror edge contributes ~ 10 to the in all cases. 



3.2.1. Power law model 

The PL fit of the spectrum (with the photon index typical 
of pulsars) suggest that the source N can be an active 
(rotation-powered) neutron star. Assuming that N is a 
pulsar^ and using the empirical relationships between the 
non-thermal X-ray and spin-down luminosities of pulsars 



(e.g. Becker & Triimper 1995, Possenti et al. 2002), one 



^ The 3.2 s individual exposure time of the ACIS instrument 
makes the search for short pulsations impossible. 



can estimate the latter, E 
is the distance to the SNR in units of 2.8 kpc. 



10^^ ^2.8 6rgs 



where d2.s 



The inferred low spin-down luminosity can be consid- 
ered as an argument against the possibility that the source 
N is a young fast-rotating pulsar with a standard mag- 
netic field (10" - 10^^ G). Therefore N can be either a 
young pulsar born with a low surface magnetic field (e.g. 
Blandford et al. 1982, Urpin et al. 1986) and/or a large 
rotation period (e.g. Spruit & Phinney 1998) or an aged 



pulsar with a decaying magnetic field. The first possibil- 
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ity would deserve detailed consideration if a candidate pe- 
riod of 12 ms found for the central X-ray source in Cas A 
(Murray et al. 2002a) is confirmed by an independent tim- 
ing analysis. The second one implies that N is an old pul- 
sar (of an arbitrary age) projected by chance on the SNR 
G 315.4-2.30. 

Note, however, that the above-mentioned empirical 
relationships are quite uncertain and should be used 
with caution. For example, the 267 ms radio pulsar PSR 
B 1853+01 associated with the SNR W44 (Wolszczan 
et al. 1991) has almost the same spin-down luminosity 
(~ 4 X lO^ergs"^) as that derived above for the source 
N, while the characteristic age of the pulsar (~ 2 x 10'* 
yr; consistent with a range of ages derived for W 4 4 usin g 
the standard Sedov- Taylor model [see Smith et al. 1985| ) 
points to its relative youth. On the other hand, the pulsar 
PSR B 1853-1-01 could be much older if its high spin-down 
rate {P ~ 0.2 x 10~*^ s s""'^) is connected to the interaction 
between the pulsar's magnetosphere and the dense circum- 
stellar mattei]^ (cf. Gvaramadze 2001 ) . The age of W 44 
also could be much greater if the origin of this mixed- 
morphology SNR is due to the cavity SN explosion (cf. 
Gvaramadze 2002). 

The uncertainties inherent to the empirical relation- 
ships are more obvious in the case of the pulsar PSR 
J 0205-^6449, recently discovered in the SNR 3C58. The 
spin-down luminosity of this young (characteristic age 
~ 5400 yr) pulsar is 2.6 x lO^^ergs""'^ (Murray et al. 
2002bD , while its (non-thermal; see Slane et al. |2002D X- 
ray luminosity is ~ 2.1 x 10'^^ c?2 e ^rgs"^, i.e. about two 
orders of magnitude less than that predicted by the em- 
pirical relationships. The disparity between the predicted 
and observed non-thermal luminosities is even higher in 
the c ase of another young pulsar - the Vela pulsar (Pavlov 
et al. 2001). Note that the photon indices and non-thermal 
luminosities of these two pulsars are almost the same as 
those predicted for the source N by the PL (or BB-I-PL) 
model, while their characteristic ages are of the same order 
of magnitude as the age of the SNR G 315.4—2.30 inferred 
in Sect. 4. Therefore one cannot exclude that the source N 
is a young pulsar with "ordinary" parameters. This should 
be tested observationally. 

3.2.2. Two-temperature blackbody model 

In the two-temperature blackbody model, the origin of the 
soft component can be attributed to the cooling surface 
of a neutron star, while the hard component to the polar 
caps heated by th e backflow of relativistic particles (e.g. 
Wang et al. [19981) . 

The best-fit BB+BB model yields the temperature of 
the soft component of 0.11 keV and the effective radius 
of 2.3 km. The inferred temperature is somewhat larger 
than that predicted by standard cooling models (e.g. Page 



1998), while the radius is significantly smaller than the 
radius of the neutron star. In principle, one can expect 
that the use of realistic neutron star atmosphere models 



(e.g. Pavlov et al. 1995) would adjust these parameters to 
acceptable values. But the interpretation of the hard com- 
ponent is more problematic. The effective polar cap area 



inferred from the best-fit BB-fBB model 



10« 



is 



too small and the observed temperature is too large to 
be consistent with models for heating of polar caps (e.g. 
Wang et al. 1998 and references therein). Therefore we 
consider this model as unplausible and suggest that at 
least at high energies the X-ray emission of the source N 
is non-thermal. 



3.2.3. Blackbody plus power law model 

The use of the BB+PL model assumes that the soft X- 
ray emission comes from the entire surface of a cooling 
neutron star or from some smaller hot areas, while the 
hard X-rays are due to the non-thermal magnetospheric 
emission. 

The best fit BB+PL model shows that the spectrum 
of the source N is dominated by the non-thermal emission 
at energies > 1 keV, and suggests the existence of a soft 
thermal component. However, the unknown interstellar 
absorption and uncertainties due to the time-dependent 
decrease in ACIS low-energy quantum efficiency make the 
estimates of the flux and bolometric luminosity of the 
BB component unconstrained. The best fit model with 
TVh as a free parameter (see Table 1) requires large in- 
terstellar absorption and consequent large bolometric lu- 
minosity. Although the large absorption, in principle, is 
consistent with values derived by Vink et al. ( 2002| ) for 
the SNR G 315.4-2.30, the large inferred bolometric lu- 
minosity implies an uncomfortably large effective radius of 
~ 60 km. The normalization of the BB component, how- 
ever, is plausible for lower values of A'^h- A more detailed 
consideration of this problem will be possible only with 
updated Chandra calibration at low energies. Note that 
dense clumps of circumstellar matter (the natural prod- 
ucts of interaction between post-main-sequence winds of 
the SN progenitor star; e.g. Gvaramadze 2001 and ref- 
erences therein; see also Sect. 4) could add a significant 
contribution to the neutral hydrogen absorption towards 
the source N, as well as possibly causing the enhanced 
absorption towards the central compact source in Cas A 
(required by some model fits of its spectrum; e.g. Murray 
et al. |2002a|). 



The large dispersion me asure variations inherent to this 
pulsar (Jacoby & Wolszczan 1999) could be considered as an 



indirect support of this suggestion. 



Note also that at present we cannot definitively con- 
firm or reject the existence of a soft BB component. 
Therefore we consider the PL model with a reasonable 
value of A^H as a good approximation for the spectrum of 
the source N. The follow-up multiwavelength observations 
of this source, including the search for pulsed emission, 
long term variability, and radio, optical or 7-ray counter- 
parts, may provide crucial information for an understand- 
ing of its nature and therefore are highly desirable. 
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4. SNR G 315.4-2.30 

We now briefly discuss a scenario for the origin of the SNR 
G 315.4—2.30 (the more detailed description will be pre- 
sented elsewhere). We note that our scenario has m uch in 
common with the model proposed by Wang et al. ( 1993 ) 
to explain the origin of large-scale structures around the 
SN 1987A, and suggest that G 315.4-2.30 is an older ver- 
sion of the latter. 

We believe that the SNR G 315.4-2.30 is the result 
of a cavity SN explosion of a moving massive star, which 
after the main-sequence (MS) phase (lasting ~ 10^ yr) 
has evolved through the red supergiant (RSG) phase 
{r-^ 10^ yr), and then experienced a short (^ 10^ yr) 
"blue loop" (i.e. the zero-age MS mass of the star was 
15 — 2OM0). During the MS phase the stellar wind (with 
the mechanical luminosity, L, of 10'^^ erg s^^) blows a 
large-scale bubble in the interstellar medium. The bubble 
eventually stalls at radius Rst ^ 12(^35/^)^/^ pc in time 
tst 10^(^35/71)1/2 where L35 = L/lO^^ergs-^ and n 
is the number density of the interstellar gas. The current 
radius of SNR G 315.4-2.30 is ~ 16 pc. The motion of that 
star causes it to cross the stalled bubble and start to in- 
teract directly with the unperturbed interstellar gas. This 
happens at time ta- ~ Rst/v* = 2.4x 10^(L35/n)i/2ti~5 yr, 
where w*^5 is the stellar velocity in units of 5kms^i. The 
time is in agreement with the duration of the MS phase if 
Vi, ~ Ikms""'^. In this case the SN progenitor star enters 
in the RSG phase while it is near the edge of the MS bub- 
ble. During the relatively short RSG phase, the star loses 
most (two thirds) of its initial mass in the form of a dense, 
slow wind. The interaction of the RSG wind with the in- 
terstellar medium results in the origin of a bow shock- 
like structure with a characteristic radius, r, determined 
by the relationship: M^sgvksgI^t^^'^ — n{2kT + mnv^) , 
where Mrsg a-nd wrsg are, correspondingly, the mass- 
loss rate and wind velocity during the RSG phase, T is 
the temperature of the ambient interstellar medium, k is 
the Boltzmann constant, and toh is the mass of a hy- 
drogen atom. For n ~ lcm~'^ (Smith 1997, Ghavamian 
" 10-5 Mq yr" 

= Ikms-i, 



et al 
T = 



2OOII ), Mrsg = 
8000 K, and 



, wrsg 
one has 



= 10kms-\ 
r ~ 1.5 pc. 



This value is in a comfortable agreement with the radius 
of the hemispherical optical nebula in the southwest of 
G 315.4-2.30 (~ 1.6 pc). 

About ^ 10'' yr before the SN explosion the progeni- 
tor star becomes a blue supergiant whose fast wind sweeps 
the material of the RSG wind. We speculate that at the 
moment of SN explosion the blue supergiant wind was 
trapped in the southwest direction by the dense material 
of the bow shock-like structure, while in the opposite di- 
rection it was able to break out in the low-density MS 
bubble, so that the SN explodes inside a "hollow" hemi- 
spherical structure open towards the MS bubble. 

The SN blast wave expands almost freely across the 
low-density MS bubble (for the explosion energy of lO^-^ 
erg and mass of the SN ejecta of ~ 3.5 M0, the expan- 
sion velocity of the blast wave is ■^bi — 5 000 km s-^) and 



reaches the northeast edge of the bubble after ~ 4.7 x 10'^ 
yr. The density jump at the edge of the bubble results in 
the abrupt deceleration of the blast wave to a velocity of 
~ {l3n\,/ny/^Viji ~ 600 — 800 kms-^ (the values derived 
from studies of Balmer-dominated filaments encircling the 
SNR G 315.4-2.30; Long fc Blair and Smith |l99^, 
where /? ~ 5 (e.g. Sgro 1975). This deceleration implies 
a density jump of a factor of ~ 200 — 300 or a number 
density of the MS bubble gas ~ 0.003 - 0.005 cm'^, 
i.e. a reasonable value given the mass lost during the MS 
phase is ~ 0.5 Mq. 

In the southwest direction, however, the expansion of 
the SN blast wave is hampered by the dense hemispheri- 
cal circumstellar shell. The shocked remains of this shell 
are now seen as the bright southwest protrusion (shown 
in Fig. ^). The existence of radiative filaments in this cor- 
ner of the SNR implies that the blast wave slows down 
to the velocity of ~ 100 kms-^, that in its turn implies 
the density contrast of ~ 15 000. For Uh — 0.003 cm"^, 
one has the number density of the circumstellar material 
of ~ 50 cm"^, i.e. in a good agreement with the density 
estimate derived by Leibowitz & Danziger ( 1983 ). On the 
other hand, the recent discovery of Balmer-dominated fil- 
aments protruding beyond the radiative arc (see Fig. 3 of 
Smith |1997| and Fig. 4 of Dickel et al. |200l| ) suggests that 
the SN blast wave has partially overrun the clumpy cir- 
cumstellar shell (cf. Franco et al. 1991) and now prop- 
agates through the interstellar medium. We believe that 
this effect is responsible for the complicated appearance 
of the southwest corner of the SNR. For the radial ex- 
tent of protrusions of c± 2 pc, and assuming that their 
mean expansion velocity is 700 km s^^, one has that 
the blowouts occured ^ 3 x 10"^ yr ago. 

5. Conclusion 

We have analyzed the archival Chandra X-Ray 
Observatory data on G 315.4—2.30 to search for a 
stellar remnant in the southwest corner of this SNR. 
The search was motivated by the hypothesis that the 
SNR G 315.4—2.30 is the result of an off-centered cavity 
SN explosion of a moving massive star, which ends its 
evolution just near the edge of the main-sequence wind- 
driven bubble. This hypothesis implies that the southwest 
protrusion in G 315.4—2.30 is the shocked material of a 
pre-existing circumstellar structure and that the actual 
location of the SN blast center is near the center of this 
structure. We have discovered two point X-ray sources 
in the "proper" place. One of the sources is interpreted 
as a foreground active star of late spectral type, while 
the second one as a candidate neutron star (perhaps a 
young "ordinary" pulsar). The follow-up observations of 
these sources will help us to understand their nature and 
thereby to test the hypothesis of the origin of the SNR 
G315.4-2.30. 

Acknowledgements. We are grateful to J. R. Dickel (the referee) 
for useful suggestions and comments. 
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